another way to view basins such as, for 
example, the yellow basin of the pendulum 
example is as if it were a lake (the basin cell) 
with three rivers draining into the lake. 
This intriguing phenomenon of Wada 
basins can be found in many applications. 
Many physical, biological, or economic sys- 
tems can be described by differential equa- 
tions analogous to those of the pendulum. 
For various choices of parameter values for 
the system, there are likely to be several 
coexisting basins. If the boundaries are 
fractal, it is likely that basin cells can be 
found. Here we showed three coexisting 
basins, but there is no limit to the number 


of coexisting basins that are Wada basins, 
all with the same boundary. Such bound- 
aries are uncertain; every boundary point is 
arbitrarily near every basin. Every boundary 
point can be perturbed arbitrarily slightly 
into any of the basins. 
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Global Patterns of Linkage 
Disequilibrium at the CD4 Locus 
and Modern Human Origins 
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Haplotypes consisting of alleles at a short tandem repeat polymorphism (STRP) and an 
Alu deletion polymorphism at the CD4 locus on chromosome 12 were analyzed in more 
than 1600 individuals sampled from 42 geographically dispersed populations (13 Af- 
rican, 2 Middle Eastern, 7 European, 9 Asian, 3 Pacific, and 8 Amerindian). Sub-Saharan 
African populations had more haplotypes and exhibited more variability in frequencies 
of haplotypes than the Northeast African or non-African populations. The Alu deletion 
was nearly always associated with a single STRP allele in non-African and Northeast 
African populations but was associated with a wide range of STRP alleles in the 
sub-Saharan African populations. This global pattern of haplotype variation and linkage 
disequilibrium Suggests a common and recent African origin for all non-African human 


populations. 





Several models for the origin of Homo sa- 
piens sapiens have been proposed. The “mul- 
tiregional origin” model suggests that there 
was no single origin for all modern humans 
(1, 2). After the radiation of Homo erectus 
from Africa into Europe and Asia 800,000 
to 1.8 million years ago (3), there was a 
continuous transition among regional pop- 
ulations from H. erectus to H. sapiens. Such 
“parallel evolution” among geographically 
dispersed populations could have been 
achieved by considerable amounts of gene 
flow between populations (1, 2). By con- 
trast, the “out of Africa” model suggests 
that all non-African human populations de- 
scend from an anatomically modern H. sa- 
piens ancestor that evolved in Africa approx- 
imately 100,000 to 200,000 years ago and 
then spread and diversified throughout the 
rest of the Earth, supplanting any Homo 
populations still present outside of Africa (1, 
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4). Migration out of Africa may have oc- 
curred in a single or in multiple waves (5). 

The best-known genetic evidence used 
to support the out of Africa hypothesis has 
come from studies of mitochondrial DNA 
(mtDNA) in which it was proposed that 
all modern mtDNA can be traced back 
through the maternal lineage to a single 
ancestor that existed in Africa between 
100,000 and 300,000 years ago (6, 7). The 
analysis and interpretation of these data 
have continued to be debated (8). Recent 
mtDNA (9) and Y chromosome (10) stud- 
ies support the original findings of a recent 
origin of all modern humans. We present 
data from the nuclear autosomal genome 
that strongly support the out of Africa 
model of human origins and provide a 
different and independent estimate, based 
on linkage disequilibrium, of the recency 
of the emigration from Africa. 
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Genetic Systems Studied 


We studied alleles from two tightly linked 
markers, located ~9.8 kb apart, within non- 
coding regions of the CD4 gene on the 
short arm of chromosome 12 (11—13) (Fig. 
1). These polymorphic markers are of two 
types that evolve with differing rates. The 
first is a short tandem repeat polymorphism 
(STRP). This class of markers consists of 
tandemly repeated blocks of two to five 
nucleotides; STRPs often have multiple al- 
leles (defined by the number of repeats) and 
moderate to high mutation rates (1/4). 
Many researchers consider them particular- 
ly useful as markers for reconstructing re- 
cent evolutionary history (15). The STRP 
at the CD4 locus consists of the pen- 
tanucleotide sequence TTTTC repeated 
between 4 and 15 times (12, 13); the prod- 
ucts (including flanking sequence) of the 
polymerase chain reaction (PCR) range in 
size from 80 base pairs (bp) for a 4-repeat 
allele to 135 bp for a 15-repeat allele (16). 
Most of the 12 alleles seen in humans are 
found primarily in Africa. Outside of Africa 
only three alleles (the 85-, 90-, and 110-bp 
alleles) ever occur at a frequency greater 
than 10%. Genotype frequencies for all 
populations are close to predicted Hardy- 
Weinberg expectations. We have also am- 
plified the CD4 STRP in common chim- 
panzee (n = 22), pygmy chimpanzee (n = 
5), gorilla (n = 5), orangutan (n = 3), and 
gibbon (n = 4). Most hominoid species are 
polymorphic, but alleles range only from 
three to six repeats (75 to 90 bp) (17). 
The second polymorphism results from 
the deletion of 256 bp of a 285-bp Alu 
element (Fig. 1) (13). This type of mutation 
is unlikely to have occurred more than 
once; DNA sequence analysis of several Alu 
deletion chromosomes from African and 
non-African individuals (18) revealed that 
all chromosomes contain the identical dele- 
tion, so that common ancestry can be as- 
sumed. The Alu deletion allele [Alu(—)| 
was typed through use of published primers 
and protocols (13); it was found to be rare or 
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absent in nine Asian, eight New World, 
and three Pacific island populations exam- 
ined and occurs at approximately a 25 to 
30% frequency in the seven European and 
two Middle Eastern populations exam- 
ined. In 13 African populations the fre- 
quency of the Alu(—) allele ranges from 7 
to 28%. Genotype frequencies for all pop- 
ulations are close to predicted Hardy- 
Weinberg expectations. The Alu(—) al- 
lele has not been detected in chimpanzees 
(n = 23), gorillas (n = 2), orangutans (n 
= 3), or gibbons (n = 4), indicating that 
the full-length Alu[Alu(+)] is the ances- 
tral state of this polymorphism and that 
the deletion event likely occurred after 
the divergence of humans from the great 
apes 4 to 6 million years ago (19). 


Haplotype Analysis 


Both markers were typed in more than 1600 
unrelated individuals from 42 geographical- 
ly dispersed populations (Table 1) (20). 
The particular combination of alleles found 
together on a single chromosome is called a 
haplotype. Haplotype frequencies were esti- 
mated either with the HAPLO program 
(21) or by simple gene counting when am- 
biguous samples were haplotyped directly 
with an allele-specific primer and long- 
range PCR (22) (Fig. 1). Mean haplotype 
frequencies and haplotype heterozygosity 
levels of populations grouped by geographic 
region are shown in Table 1. Haplotype 
frequencies for: several individual popula- 
tions are represented graphically in Fig. 2. 
Only the 85-bp Alu(+) and the 110-bp 
Alu(+) haplotypes occur at a frequency 
greater than 10% in all non-African popu- 
lations. The 90-bp Alu(—) haplotype oc- 
curs at a high frequency (23 to 34%) in the 
Middle Eastern and European populations 
and is rare or absent in all Asian, Pacific 
island, and New World populations (0 to 
9% frequency) (Fig. 2, A to D). Among 


non-African populations heterozygosity is 
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highest in the Middle East (0.75 in the . World (0.37 in the Brazilian Ticuna) and 
Israeli Druze), decreases from west to east Australo-Melanesia (0.28 in the Nasioi 
across Eurasia, and is lowest in the New Melanesians) (Table 1). 


Fig. 1. Location and dis- 
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show the location and di- 
rection of primers used for amplification by PCR of the STRP (A and B) (72) and Alu polymorphism (C and 
D) (73). Primer E (5’-GCGCAAGATGCACAGCCTG-3’) is designed to be specific for the Alu(—) allele such 
that long-range PCR amplification with primer A occurs only on chromosomes that contain the Alu 
deletion (22). A second round of amplification with primers A and B with 1 jl of the 9.9-kb PCR product 
allows determination of the size of the STRP allele on the Alu(—) chromosome. 
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Fig. 2. CD4 STRP-Alu(+/—) haplotype frequencies. Solid blue bars, Alu(+) chromosomes; striped red 
bars, Alu(—) chromosomes. Haplotype frequencies are shown for four diverse non-African (A to D) and 
four diverse African (E to H) populations selected from the geographic regions shown in Table 1. 
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Haplotype frequencies in Egyptians do 
not differ significantly from those in Mid- 
dle Eastern populations (23), consistent 
with linguistic evidence suggesting a re- 
cent common ancestry. The Ethiopian 
Jewish (Fig. 2E) and Somali populations 
show a pattern of haplotype frequencies 
similar to, but significantly different from, 
the non-African populations, with slightly 


Table 1. Unweighted mean CD4 STRP-Alu frequency estimates and standard 
errors (in parentheses) for seven regions of the world. Haplotypes are represented 
by the STRP allele and the Alu allele. 2, number of populations; 2N, number of 


more haplotypes (mean heterozygosity of 
0.80) (24). 

The pattern of haplotype frequencies in 
sub-Saharan African populations differs from 
that in the Northeast African and non-Af- 
rican populations (Fig. 2, F to H). Compared 
to non-African populations, the sub-Saharan 
populations have more haplotypes that exist 
at a higher frequency; nine haplotypes occur 


STRP (bp)-Alu(-+) haplotypes 


at a frequency greater than 10% in at least 
one population, and heterozygosity ranges 
from 0.82 to 0.90. Frequencies of haplotypes 
vary considerably among populations as il- 
lustrated by the Kikuyu, Woloff, and Nama 
populations (Fig. 2, F to H). 

The nonrandom association between 


the CD4 STRP alleles and the Alu(—) 


allele is referred to as linkage disequilibrium 


chromosomes sampled. No 135-bp Alu(—) haplotypes were observed. Estimated 
heterozygosities (Het.) for the regions were calculated for each population as 1 - 
>[(frequency of each haplotype)*] and then averaged across populations. 


n ON 80+ 85+ 90+ 95+ 100+ 105+ 110+ 115+ 120+ 125+ 130+ 135+ 
New World* 
8 (04 O Qs553 0.001 O 0.001 0.001 0.421 0.004 0.001 O O O 
(0.047) (0.001) (0.001) (0.001) (0.051) (0.004) (0.001) 
Pacific island and Australo-Melanesiant 
e) 130 O 0.727 0.013 O O 0.007 O27 0.013 O O O O 
(0.094) (0.007) (0.007) (0.088) (0.007) 
Asiant 
9 600 O 0.603 0.011 0.003 0.003 0.004 O30 0.023 0.002 O 0 O 
(0.038) (0.005) (0.002) (0.003) (0.003) (0.033) (0.011) (0.001) 
Europeans 
4 658 O 0.361 0.010 0.003 O 0.007 0.289 0.037 0.017 O O O 
(0.038) (0.003) (0,002) (0.006) (0.030) (0.009) (0.005) 
Middle Eastern|| 
2 198 0 0.335 0.015 0 O O 0.240 0.045 0.020 0.025 O O 
(0.075) (0.005) (0.030) (0.005) (0.000) (0.025) 
Northeast African 
3 232 O O23 0.010 O O 0.020 0.280 0.113 0.063 0.006 0.027 O 
(0.058) (0.010) (0.012) (0.026) (0.047) (0.038) (0.006) (0.007) 
Sub-Saharan African# | 
10 806 0,002 0.187 0.023 0.030 0.074 0.026 0.172 0.156 0.107 0.024 0.021 0.005 
(0.002) (0.022) (0.012) (0.013) (0.01 7) (0.009) (0.026) (0.026) (0.019) (0.007) (0.007) 0.002) 
STRP (bp)—Alu(—) haplotypes 
n ON 80— 85— 90 — 95— 100— 105-— 110- 115- 120- 125- 130— Het. 
New World” 
8 (04 O O 0.018 0 O O O O O 0 O 0.48 
(0.010) (0.018) 
Pacific island and Australo-Melanesiant 
S) 130 O O 0.007 O O O O O O O O 0.40 
(0.007) (0.080) 
Asiant 
9 600 O 0.001 0.037 0 O O O O O O O 0.52 
(0.001) (0.012) (0.034) 
Europeans 
i 658 O 0.001 0.276 0.001 O O 0.001 O O O O 0.69 
(0.001) (0.017) (0.001) (0.001) (0.016) 
Middle Eastern|| 
2 198 O 0.005 0.305 O 0.005 O O O O O O 0.73 
(0.005) (0.015) (0.005) (0.025) 
Northeast African 
3 232 O 0.020 0.160 0.010 0.004 0.013 O O O O O 0.78 
(0.012) (0.044) (0.010) (0.008) (0.013 (0.025) 
Sub-Saharan African# 
10 806 O 0.046 0.042 0.009 0.007 0.003 O 0.054 0.007 0.005 0.001 0.85 
(0.011) (0.008) (0.003) (0.003) (0.003) (0.014) (0.004) (0.002) (0.001) (0.010) 


*Ticuna (Brazil; 2N = 124), Karitiana (Brazil; 2N = 100), Surui (Brazil; 2N = 78), Quechua (Peru; 2N = 44), Maya (Yucatan; 2N = 98), Jemez Pueblo (USA; 2N = 86), Pima (USA; 2N 
= 70), Cheyenne (USA; 2N = 104). tMicronesian (2N = 46), New Guinean (2N = 40), Nasioi (Bougainville; 2N = 44). Malaysian (2N = 24), Atayal (Taiwan; 2N = 74), Ami 
(Taiwan; 2N = 68), Kochari (Assam; 2N = 32), Cambodian (2N = 48), Chinese (2N = 88), Japanese (2N = 86), Yakut (Siberia; 2N = 72), Asiatic Indian (2N = 108). §Adygei 
(Caucasus; 2N = 98), Europeans (mixed; 2N = 178), Finns (2N = 68), Danes (2N = 38), Basque (Spain; 2N = 124), Ashkenazi Jews (2N = 102), Roman Jews (2N = 50). Druze 
(Israel; 2N = 96), Yemenite Jews (2N = 102). qEgyptian (2N = 68), Ethiopian Jews (2N = 114), Somali (2N = 50). #Kikuyu (Kenya; 2N = 44), Yoruba (Nigeria; 2N = 36), Wolof 
(Senegal; 2N = 102), Mbuti (Zaire; 2N = 74), Biaka (Central African Republic; 2N = 106), Herero Bantu-speakers (Namibia; 2N = 98), Zu/Wasi !Kung San (Namibia; 2NV = 88), Sekele 
[Kung San (Namibia; 2N = 104) Bantu-speakers (mixed South African; 2N = 88), Nama (Namibia; 2N = 66). 
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and can be estimated from the parameter 8 
(25), where 8(_.) = [pc_)(x) — Pea yOI/1 
— poay(x)], pr_)(x) is the proportion of 
STRP allele x on Alu(—) chromosomes, 
and p,, (x) is the proportion of STRP allele 
x on Alu(+) chromosomes (Table 2). Out- 
side of Africa only one value is large: 8,_ 99) 
= 0.98 for combined non-African popula- 
tions, a value close to the maximum of 1.0, 
indicating numerically the obvious finding 
(Table 1) that the Alu deletion allele al- 
most always occurs on a chromosome with 
the 90-bp STRP allele. The largest 8,_ 
values for the Egyptian, Ethiopian, and So- 
mali populations are also with the 90-bp 
STRP allele [8,_ 99) values of 1.0, 0.61, and 
0.67, respectively). Among the sub-Saharan 
African populations, examples can be found 
of the Alu(—) allele positively associated 
with several STRP alleles. However, only 
the 90-bp STRP allele consistently (in 9 


outof 10 populations) exhibits a positive 
association with the Alu(—) allele; it also 
has the highest median 8,_ ,) value (0.24), 
suggesting that the Alu Jeetion originally 
occurred on a chromosome containing a 
90-bp STRP allele and there has not been 
sufficient time for complete equilibrium to 
be attained. Figure 3 summarizes the three 
categories of 8(_ 99) values: strong in non- 
Africans (and Egyptians), intermediate in 
Ethiopian and Somali populations, and still 
present but weak in sub-Saharan African 
populations. 


Table 2. 5._ 49) values for the three CD4 STRP 
alleles most frequently associated with the Alu(—) 
allele. African populations are listed individually, 
and all non-African populations containing the 
Alu(—) allele are combined for a single estimate. A 
value of +1 indicates complete association with 
the Alu(—) allele, O indicates no association, anda 
dash indicates that the STRP allele is more com- 
mon on Alu(+) than on Alu(—) chromosomes. No 
other STRP allele had a 8_ ,, value greater than 
0.20 in any population. 


STRP allele (bp) 


Population —————— 
85 90 115 
Non-African 
Combined - 0.98 - 
Northeast African 

Egyptian — 1 — 

Ethiopian Jewish — 0.61 — 

Somalian 0.18 0.67 - 

Sub-Saharan African 

Woloff 0.41 0.22 _ 
Yoruba - 0.20 0.18 
Mbuti 0.46 0.19 0.00 

Biaka. 0.38 - - 
Kikuyu - 0.11 0.64 
Bantu-speakers - 0.26 0.29 
Herero 0.01 0.38 0.13 
Nama - 0.27 0.18 
Sekele !Kung San - 0.26 0.13 
Zu/wasi !Kung San - 0.19 0.03 


Implications of Linkage 
Disequilibrium 


Distributions of haplotype frequencies and 
linkage disequilibrium in extant popula- 
tions are the result of several processes: 
mutation at the STRP, recombination be- 
tween the two sites, random genetic drift, 
and gene flow among populations. We be- 
gin our analysis by reconstructing a plausi- 
ble scenario to explain the contrasting pat- 
terns of linkage disequilibrium in Africans 
and non-Africans. 

The CD4 STRP apparently started as a 
small-sized repeat, still seen in nonhuman 
primates, and has accumulated variation in 
the human lineage during the 4 to 6 million 
years since humans diverged from the great 
apes (19). At some time during that period 
the Alu was partially deleted from a single 
chromosome which, through genetic drift, 
eventually became common. Initially, there 
must have been complete linkage disequi- 
librium between the Alu(—) allele and the 
“ancestral” nearby STRP allele (most likely 
a 90-bp allele, as noted above). Over time, 
because of the effects of mutation and re- 
combination, that linkage disequilibrium 
was largely lost in the ancestral African 
population. 

Then, through random genetic drift, the 
Alu(—) chromosome containing the 90-bp 
STRP allele became the predominant 
Alu(—) chromosome in one or more popu- 
lations in the Northeastern periphery of 
Africa. Further reduction in genetic varia- 
tion most likely occurred when members of 
this ancestral population left Africa and 
populated the rest of the world, bringing 
with them the chromosome containing the 
90-bp STRP allele and the Alu(—) allele. 
The presence of almost complete linkage 
disequilibrium between these two alleles in 
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all non-African populations surveyed sug- 
gests that the Alu(—) chromosome con- 
taining the 90-bp STRP allele was the only 
Alu(—) chromosome to remain in the pop- 
ulations that left Africa. The migration from 
Africa was also accompanied by a further 
reduction in the number of Alu(+) haplo- 
types; only those Alu( +) chromosomes con- 
taining an 85- or 110-bp STRP allele remain 
at a significant frequency outside of Africa. 
This finding supports studies of both 
mtDNA and nuclear DNA, suggesting that 
the ancestral population that left Africa was 
small (26, 27), but contrasts with results 
from the major histocompatibility complex 
(MHC) locus, for which there remains con- 
siderable variation outside of Africa (28). 
These contrasting patterns may be due to 
selection at the MHC locus that favors 
maintenance of allelic lineages (29). 

The intermediate pattern of genetic 
variation seen in the Ethiopian Jewish and 
Somali populations could be explained ei- 
ther by recent admixture with nearby Mid- 
dle Eastern populations or a recent com- 
mon ancestry with non-African popula- 
tions (30). Because these two Northeast- 
ern African populations have only a subset 
of the haplotype diversity present in sub- 
Saharan Africa, admixture between sub- 
Saharan and Middle Eastern populations 
does not seem a good explanation of their 
haplotype frequencies; rather, these popu- 
lations may represent the modern descen- 
dants of the ancestral population that 
spawned the migration from Africa. 

The high level of haplotype diversity 
observed within all 10 sub-Saharan African 
populations of diverse geographic origin 
suggests that much of the variation ob- 
served today was present in the ancestral 
African population from which they de- 
scended, with differences among these pop- 





Fig. 3. Global distribution of populations used in the CD4 STRP-Alu haplotype study and the associated 


3, 90) Values for the 90-bp STRP and Alu(—) alleles. (®) Populations with no Alu(—) allele detected; (A) 
populations with 8,_ 99) values ranging from 0.90 to 1.0; (Ml) populations with 8/_ gq) values ranging from 
0.63 to 0.80; (@) populations with 8._ 99) values less than 0.38. 
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ulations due to recent effects of random 
genetic drift in each region. Thus, these 
data would suggest that sub-Saharan Afri- 
cans have maintained a large effective pop- 
ulation size and that no recent bottleneck 
events have occurred during the evolution 
of H. sapiens in Africa. 

The much greater variation in both 
number and frequency of CD4 STRP-Alu 
haplotypes in sub-Saharan Africans is con- 
sistent with the out of Africa hypothesis 
(4). The clinal pattern of decreasing het- 
erozygosities west to east outside of Africa is 
consistent with a migration event of mod- 
ern humans out of Northeast Africa into 
the Middle East and Europe, and east into 
Asia, the Pacific islands, and the New 
World. In addition, all non-African popu- 
lations share a similar pattern of haplotype 
variation and linkage disequilibrium, which 
implies that they all derive from a single 
ancestral gene pool. This pattern is similar 
to that seen in Northeast African popula- 
tions but is distinct from that in all sub- 
Saharan African populations. Patterns of 
genetic variation distinguishing between 
African and non-African populations have 
been documented in data of both mtDNA 
(6, 31) and nuclear DNA (27, 32, 33). In 
addition, however, the retention of essen- 
tially complete linkage disequilibrium out- 
side of Africa implies that the common 
origin of the non-African populations is 
recent because there has not been sufficient 
time for linkage disequilibrium to be dimin- 
ished by the shuffling effects of mutation 
and recombination. 





Dating the Emergence 
from Africa 


Assuming the above scenario, we can at- 
tempt to estimate how long ago the ances- 
tors of modern non-Africans existed as a 
single distinct population with complete 





Table 3. Relative weighted mean frequencies of the CD4 STRP alleles in 
each region standardized separately for Alu(+) and Alu(—) chromosomes. 


linkage disequilibrium of the 90-bp and 


Alu(—) alleles. All Alu(—) alleles appear / 





to be descended from a single original mu- 
tation event that occurred in the human 
lineage a maximum of 4 to 6 million years 
ago (19). There has been enough time in 
sub-Saharan Africa since the origin of the 
Alu deletion for mutation and recombina- 
tion to increase the allelic variation at the 
STRP on Alu(—) chromosomes, but not as 
much time since the introduction of the 
Alu(—) chromosome outside of Africa for 
variation to be regenerated. 

To assess the relative effect of recombi- 
nation and mutation on reducing levels of 
linkage disequilibrium between the Alu de- 
letion and the STRP, we contrast the dis- 
tributions of STRP alleles on Alu(+) and 
Alu(—) chromosomes in sub-Saharan Afri- 
can populations. High levels of recombina- 
tion will make those distributions similar 
within populations. For comparison, we 
contrast the distributions of STRP alleles 
on chromosomes with the same Alu allele 
from different sub-Saharan populations. 
Genetic drift will tend to make those dis- 
tributions different, whereas gene flow will 
tend to make them similar. By any of sev- 
eral statistical methods (34) we see high 
levels of similarity among Alu(+) chromo- 
somes and among Alu(—) chromosomes 
from different populations and much less 
similarity between Alu(+) and Alu(—) 
chromosomes taken from the same or from 
different populations. These results imply a 
recent divergence of the populations or 
considerable gene flow among them (or 
both) and suggest that recombination is 
playing only a negligible role in equalizing 
the allele frequency distributions on 
Alu(+) versus Alu(—) chromosomes with- 
in populations. This conclusion is supported 
by the observation that no Alu(—) chro- 
mosomes (out of 132) in Sub-Saharan Af- 
rican population and only 1 chromosome 

















2N, number of chromosomes sampled. Haplotype diversity d(x) was calcu- 


ON 80+ 85+ 90+ 95+ 
2020 «OO 0.558 0.010 0.002 
674 0.001 0.203 0.028 0.040 

ON 80— 85— o0= 95— 
270 0 0.011 0.978 0.004 
132 0 0.264 0.258 0.068 
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STRP (op)—Alu(+) haplotypes 


(out of 286) outside of Africa carries a 
110-bp allele, despite this allele having a 
frequency of 21% on Alu(+) chromosomes 
across Sub-Saharan Africa and 39% outside 
of Africa (Table 3). In addition, most of the 
variance of STRP alleles on Alu(—) chro- 
mosomes in the Northeast African and 
non-African populations consists of alleles 
smaller than 110 bp. This pattern of varia- 
tion is consistent with a model of mutation 
involving primarily small changes in repeat 
size away from the ancestral 90-bp allele 
(35). Thus, both inside and outside Africa, 
the level of linkage disequilibrium now ob- 
served appears to result primarily from mu- 
tation occurring over time, with recombi- 
nation having a smaller effect. 

Random fluctuations in allele frequen- 
cies—random genetic drift—confound ef- 
forts to make precise estimates of how many 
generations have elapsed since complete 
disequilibrium existed. We cannot estimate 
precisely when the founding event out of 
Africa occurred without knowing the mu- 
tation and recombination rates, which are 
too small to be estimated by direct observa- 
tion. However, by comparing the relative 
breakdown of disequilibrium in sub-Sahar- 
an and non-African populations, we can 
make a conservative estimate of the maxi- 
mum time since the last common ancestor 
of the non-African populations (36). For 
our analysis we make the following assump- 
tions: (1) The Alu deletion originally oc- 
curred in Africa as a single event after the 
separation of humans from the great apes 
and was originally associated with a single 
STRP allele (the 90-bp allele). (2) In the 
founding population for non-Africans, the 
Alu(—) chromosome was associated with 
only a single STRP allele (the 90-bp allele). 
(3) The mutation rate at the STRP is the 
same inside and outside of Africa. (4) There 
has been no back mutation to a 90-bp allele 
after a mutation away from the 90-bp allele. 














lated separately for Alu(+) and Alu(—) chromosomes as 1 — f* where fis the 
standardized frequency given in the table for STRP allele /. 


100+ 105+ 110+ 115+ 120+ 125+ 130+ 1385+ a[Alu(+)] 
Non-African 
0.002 0.005 0.390 O;027 0.009 0.002 O O Oro4 
Sub-Saharan African 
0.090 0.030 0.210 0.200 0.140 0.0380 0.030 0.004 0.84 
STRP (op)—Alu(—) haplotypes 
100-— 105- 10S a 120-— 125- 130— 135- d[Alu(—)] 
Non-African 
0.004 O 0.004 O O O O O 0.04 
Sub-Saharan African 
0.041 0.008 O 0271 0.045 0.0388 0.007 O 0.78 
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(5) Except for the founding “event,” muta- 
tion has had a larger effect than genetic 
drift or recombination on the proportions of 
most STRP alleles on Alu(+) and Alu(—) 
chromosomes (37). 

Let P, represent the frequency of the 
progenitor STRP allele on non-African 
Alu(—) chromosomes. Assuming a muta- 
tion rate of w per generation, the frequency 
of chromosomes still bearing the original 
mutation can be approximated by e—®#, 
where Nz, is the number of generations 
since the founding. Similarly, within Afri- 
ca, if P, is the frequency of chromosomes 
still bearing the original STRP allele, P, 
can be approximated by e—N4". Thus, to 
obtain the relative age R of the Alu(—) 
chromosome inside versus outside of Africa, 
we can calculate 


_ In (Pa) _ Na 
~ in (Py) Ny o 





The relative frequencies of STRP alleles 
separately on Alu(+) and on Alu(—) chro- 
mosomes are given in Table 3 for the 
weighted averages of sub-Saharan and of 
non-African populations (38). Using the 
pooled frequencies in each region minimiz- 
es the effects of recent population-specific 
genetic drift. Outside Africa it is clear that 
the 90-bp allele was the progenitor on 
Alu(—) chromosomes, and its current fre- 
quency on such chromosomes is 264/270 = 
0.978. Within sub-Saharan Africa, three 
alleles occur at nearly equal high frequen- 
cies on Alu(—) chromosomes—the 85-, 
90-, and 115-bp alleles. As previously not- 
ed, patterns of linkage disequilibrium sug- 
gest that the ancestral Alu(—) chromosome 
in Africa contained a 90-bp STRP allele 
(Table 2). The high frequency of the 85-bp 
allele on Alu(—) chromosomes may be the 
result of genetic drift or of selective or 
mechanistic constraint against mutation to 
alleles smaller than 85 bp. Sequencing anal- 
ysis of several 115-bp alleles from Alu(—) 
chromosomes suggests that their presence 
on Alu(—) chromosomes results from an 
ancient recombination event with an 
Alu(+) chromosome and _ they have 
reached high frequency through genetic 
drift (39). Because we restrict our analyses 
to include only mutation events from the 
ancestral 90-bp allele and to make our anal- 
yses as conservative as possible, we exclude 
from initial consideration all STRP alleles 
of 110 bp or greater (1 case outside Africa, 
49 cases within Africa). We find P, = 
34/85 = 0.40 and P, = 264/269 = 0.9814. 
Therefore, R = In(P,)/In(Pg) = 48.8. As- 
suming a maximum of 5 million years before 
present (Y.B.P.) for the origin of the 
Alu(—) chromosome in Africa, this gives a 
maximum age of 102,000 Y.B.P. [with an 
upper bound of 217,000 to 313,000 Y.B.P. 








(40)] for the introduction of the Alu(—) 
chromosome outside of Africa. Slightly dif- 
ferent approaches to analyzing the data give 
somewhat -different values of R and hence 
different dates for the founding of the non- 
African populations (41). However, even 
with implausibly conservative assumptions, 
the dates that result from our calculations are 
much more recent than the dates proposed 
by the multiregional origin hypothesis. 

Our date estimates represent maximum 
values based on assumption | enumerated 
above, that the Alu(—) allele arose 5 million 
years ago. It is likely that the Alu deletion 
event occurred more recently, in which case 
our estimates of the date of founding of the 
non-African populations would also be more 
recent (42). Assumption 2 also gives a con- 
servative estimate because violation of this 
assumption would result in a larger value of 
R. Regarding assumption 4, back mutation 
outside of Africa is unlikely to have played a 
large role because of the very low frequency 
of alleles other than 90 bp. Back mutation 
having occurred within Africa is more likely; 
however, ignoring such back mutation leads 
to an underestimate of N,, and thus R. 
Hence, assumption 4 also provides a conser- 
vative estimate of R. 


A Single Recent African Origin 


Because this is a single locus subject to the 
effects of random genetic drift and we are 
observing only one realization of evolution, 
it is likely that similar analysis of other loci 
will give slightly different date estimates. 
However, the difference in haplotype vari- 
ation at the CD4 locus between sub-Sahar- 
an African and non-African populations is 
striking and cannot easily be accounted for 
except by a recent common origin of non- 
Africans from Africa. We present conserva- 
tive maximum estimates for the origin of 
non-African humans that are consistent 
with the fossil record, which places anatom- 
ically modern humans in Africa as early as 
120,000 Y.B.P. and dates the earliest modern 
human fossils in the Middle East at 90,000 to 
120,000 Y.B.P. (43). They are also consist- 
ent with the estimates for the date for the 
last common mtDNA ancestor of all modern 
humans (African and non-African) of 
100,000 to 300,000 Y.B.P. (6, 9). 

Two alternative hypotheses that might 
explain the global pattern of haplotype 
variation and linkage disequilibrium are se- 
lection and the multiregional origin hy- 
pothesis. Specific selection schemes acting 
on mutations closely linked to the CD4 
STRP can either reduce the variance of 
alleles at the STRP (a selective “sweep”) or 
result in higher than expected levels of 
variation at the STRP (overdominant se- 
lection) (44). However, no simple selective 
schema can explain the diversity of haplo- 
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types in sub-Saharan Africa and the specific 
pattern of haplotypes seen in all non-Afri- 
can populations (45). Similar patterns of 
haplotype variation and differential linkage 
disequilibrium between Africans and non- 
Africans at several additional unlinked loci 
(33, 46) also suggest that selection effects 
could not explain these results. 

The second alternative, the multire- 
gional origin hypothesis, seems highly un- 
likely: It predicts roughly equivalent time 
depth and genetic diversity in all parts of 
the world. There has been enough time 
since the origin of the Alu deletion in 
Africa for mutation at the STRP, and to a 
lesser extent recombination between the 
markers, to generate multiple alleles and 
nearly randomize the Alu(—) allele with 
respect to the specific STRP allele on the 
same chromosome. In contrast, the limited 
number of STRP alleles and the essentially 
complete linkage disequilibrium in all non- 
African populations suggest a recent single 
population for their origin (47). It is highly 
unlikely that this pattern could be ex- 
plained by a recent bottleneck event in 
non-African populations that have been 
diverging since H. erectus emerged from 
Africa 800,000 or more years ago; we would 
not expect to find such similar patterns of 
haplotype diversity and linkage disequilib- 
rium. An out of Africa model invoking 
expansion and spreading of populations 
from a single common ancestor is most 
parsimonious. 

Because we can trace the inheritance of 
chromosomes through both males and fe- 
males, analysis of nuclear DNA markers 
provide a useful comparison with studies of 
mitochondrial and Y chromosome DNA 
variation. Haplotyped loci are especially 
useful because linkage disequilibrium pro- 
vides more information beyond that of the 
single polymorphisms considered as inde- 
pendent markers. By choosing haplotype 
systems that approach equilibrium with dif- 
ferent rates, we can explore different evolu- 
tionary time depths. The haplotype system 
at the CD4 locus has the appropriate low 
levels of recombination and mutation to be 
particularly useful for studying the period 
when modern humans first expanded out of 
Africa. The decreased variation and strong 
disequilibrium seen in non-African popula- 
tions contrast with the quite different pat- 
tern in sub-Saharan African populations 
and together constitute some of the stron- 
gest evidence to date from the nuclear ge- 
nome for a recent common African origin of 
all modern non-African human populations. 
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before rapid population expansion, we could be esti- 
mating the time since divergence and expansion of 
the ancestral non-African population, which may have 
occurred subsequent to migration out of Africa. 
Other assumptions may be warranted for other hap- 
lotype loci. Specifically, for markers with low muta- 
tion rate but proportionately higher recombination 
rates than evident at CD4, models and analyses 
based primarily on recombination would be most 
appropriate. 
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excluded from the analysis because of the possibility 
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Sequencing analysis of more than 100 alleles from 
chromosomes of humans of diverse geographical 
origin shows that alleles of 110 bp or larger on Alu(+) 
chromosomes have an imperfect CTTTC as the 
fourth repeat from the 5’ end. Thus, the ancestral 
STRP allele was most likely a small-sized perfect 
repeat of TT TTC (77), and alleles 110 bp or larger 
are derived from a single mutation that generated a 
higher repeat-number allele with the imperfect 
CT TTC. Sequencing analysis of several large-sized 
alleles from Alu(—) chromosomes demonstrates that 
they share the point mutation present on Alu(+) 
chromosomes; thus, these alleles most likely derive 
from an ancient recombination event transferring a 
large-sized allele from an Alu(+) chromosome to an 
Alu(—) chromosome with subsequent mutation by 
one or few repeats to other large-sized alleles. 

To estimate a lower bound for AR from these data 
(hence an upper bound maximum age), we note that 
the variability in the estimate of R depends much 
more on variability in the denominator than the nu- 
merator. According to a Poisson distribution, the 
value of P,, that corresponds to a probability of 5% 
for there to appear at most five alleles outside of 
Africa without the 90 bp STRP on Alu(—) chromo- 
somes is 0.961. This gives a lower 95% confidence 
limit for R of In(O.40)/In(0.961) = 23.0 (or an age 
maximum of 5 million Y.B.P./23.0 = 217,000 
Y.B.P.). The Poisson assumption for non-Africans 
depends on statistical independence of the chromo- 
somes sampled, a reasonable assumption because 
at least five of the six non- 90-bp Alu(—) haplotypes 
derive from distinct populations (one 100-bp in the 
Druze, one 85-bp in Yemenite Jews, one 85-bp and 
one 95-bp in mixed European, one 85-bp in Asiatic 
Indian, and one 110-bp in Adygei). Regarding vari- 
ability on the value of P,, we examined the P values 
for the five sub-Saharan populations with more than 
10 Alu(—) chromosomes. The mean of these values 
is 0.40, the same as the combined value we used 
above. However, even if we use the largest of these 
values (0.53), along with the lower 95% confidence 
limit for P,, we obtain a value for FA of In(Q.53)/ 
In(0.961) = 16.0 (or maximum age of 5 million/16.0 
= 313,000 Y.B.P.). 

If the 85-bp allele rose to high frequency in sub- 
Saharan Africa through drift or through mechanistic 
or selective constraint, we can combine the 85- and 
90-bp alleles into a single progenitor group to obtain 
a highly conservative estimate of A. In this case, 
using alleles less than 110 bp in size, P, = 69/85 = 
0.812; applying Eq. 1 gives R = 11.1 for a maximal 
age of 450,000 years for the origin of the Alu(—) 
chromosome outside of Africa. Assuming the 115- 
bp allele on Alu(—) chromosomes is the result of an 
early recombination event (39), the 9O- and 115-bp 
alleles can be considered as a single progenitor 
group, and STRP alleles of all other sizes as the 
mutants, giving P, = 70/132, P, = 264/270, R = 
28.2, for a maximal age of 177,000 years. Various 


models that assume stepwise mutation of the STRP 
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early with time (48) [D. B. Goldstein, A. Ruiz Linares, 
L. L. Cavalli-Sforza, M. W. Feldman, Genetics 139, 
463 (1995)]. Thus, we can estimate R as the ratio of 
variances of allele sizes (in base pairs) for alleles less 
than 110 bp (22.5/0.75 = 30.0), giving a maximum 
age of 167,000 Y.B.P. Because we have considered 
only a single locus, the standard error on this vari- 
ance-ratio estimate is high. However, the apparent 
lower boundary of 85 bp to the STRP allele size 
results in an underestimate of the time of origin of the 
Alu(—) chromosome in Africa, so again this estimate 
of FR is conservative. 
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some, it is still quite ancient. 
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“non- 90 bp”’ Alu(—) chromosomes by recombina- 
tion between the markers or mutation at the STRP. In 
addition, the very low frequency of the Alu(—) chro- 
mosomes in Asia, the Pacific islands, and the New 
World, as well as the very high frequency of 85-bp 
Alu(+) and 110-bp Alu(+) chromosomes in all non- 
African populations, argues against strong positive 
selection for the Alu(—) chromosome. 
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The Exchange of Impact Ejecta 
Between Terrestrial Planets 


Brett J. Gladman, Joseph A. Burns, Martin Duncan, Pascal Lee, 
Harold F. Levison 


Orbital histories of ejecta from the terrestrial planets were numerically integrated to 
study their transfer to Earth. The properties of the lunar and martian meteorites are 
consistent with a recurrent ejection of small meteoroids as a result of impacts on their 
parent bodies. Long-range gravitational effects, especially secular resonances, strong- 
ly influence the orbits of many meteoroids, increasing their collision rates with other 
planets and the sun. These effects and collisional destruction in the asteroid belt result 
in shortened time scales and higher fluxes than previously believed, especially for 
martian meteorites. A small flux of mercurian ejecta appears possible; recovery of 
meteorites from the Earth and Venus is less likely. 





The study of meteorites has illuminated 
the nature of extraterrestrial environments 
and astrophysical processes, particularly the 
conditions at the time of our solar system’s 
formation. Most meteorites come from as- 
teroids, but recently a number of objects 
from the moon and Mars have been recog- 
nized. These latter meteorites help to char- 
acterize the surfaces of these bodies, espe- 
cially the martian meteorites, which are our 
only samples of that planet. To learn more 
about the parent bodies and the paths that 
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the meteorites traveled before arriving on 
Earth, we must understand the orbital dy- 
namics governing their transfer. In particu- 
lar, what is the delivery efficiency, that is, 
the fraction of escaping ejecta that reach 
Earth, from different sources? 

The SNC meteorites (1) have features 
that suggest they are derived from Mars (2) 
and are somehow delivered to Earth. Even 
though the petrology and young crystalliza- 
tion ages of the SNCs point to an origin on 
a large parent body with recent geologic 
activity, Mars was only recently accepted as 
their source because it was thought unlikely 
that rocks could survive being blasted off of 
a planet (3). Another argument against a 
martian origin was that if these objects were 
launched from Mars, surely there should be 
many more meteorites coming from the 
moon, which is closer and has a lower escape 
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bp Alu(—) haplotype to have been recently intro- 
duced into preexisting H. erectus populations and 
achieve the frequencies seen today. The low fre- 
quency of the 90-bp Alu(—) chromosome in all Asian, 
Pacific island, and New World populations argues 
against high levels of gene flow from European or 
Middle Eastern populations into these regions before 
historical times. 
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velocity; yet, as of 1982, there were six SNC 
falls but no lunar meteorites. Then the lunar 
meteorite ALHA81005 was recognized 
among Antarctic meteorites (4), and be- 
cause of our familiarity with the returned 
Apollo and Luna samples, the origin of this 
meteorite was immediately accepted. With 
more Antarctic meteorites (5) being recov- 
ered, we currently have a dozen members of 
each class, although they comprise only 
0.1% of all meteorites. The SNCs are now 
generally acknowledged to have originated 
on Mars, after further examination of their 
composition, especially the virtually perfect 
isotopic match between gases trapped within 
one of them and the martian atmosphere, as 
determined by the Viking landers (6). Thus, 
we use the term “martian meteorites” for the 
SNCs and ALH84001, the latter being dis- 
tinct from the SNC classification (7). 

We can learn about the dynamics of the 
inner solar system by comparing the mea- 
sured transfer ages of such meteorites against 
orbital histories. Some of the ideas presented 
below are not new, but our dynamical sim- 
ulations improve on previous work that, be- 
cause of computational limitations of the 
time, used Monte Carlo calculations rather 
than full orbital integrations. In principle, 
some aspects of the cratering process itself 
can also be constrained (8). 

A cosmic ray exposure (CRE) age (9) of 
a meteorite is the time during which the 
object was bombarded by energetic cosmic 
rays in space. During this exposure, measur- 
able radioactive isotopes accumulate, allow- 
ing the duration of the exposure to be esti- 
mated (Table 1). Most lunar meteorites 
were delivered to the Earth in a far shorter 
time than any martian meteorite, and the 
martian meteorites have an average mass 38 
times that of the lunar ones. The issue of 
pairing (10) is unlikely to affect these 
trends significantly. 
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